Stability and instability in parametric resonance and quantum Zeno effect 
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A quantum mechanical version of a classical inverted pen- 
dulum is analyzed. The stabilization of the classical motion is 
reflected in the bounded evolution of the quantum mechani- 
cal operators in the Heisenberg picture. Interesting links with 
the quantum Zeno effect are discussed. 

PACS numbers: 42.65. Sf ; Q3.65.Bz; 05.45.-a; 42.65.Yj 

An inverted pendulum is an ordinary classical pen- 
dulum initially prepared in the vertical upright position 
10-^. This is normally an unstable system, but can be 
made stable by imposing a vertical oscillatory motion to 
the pivot. In a few words, when the pivot is acceler- 
ated upwards the motion is unstable, while when it is 
accelerated downwards the motion can be stable: the pe- 
riodic switch between these two situations can be glob- 
ally stable or unstable depending on the values of some 
physical parameters. In particular, when the frequency 
of the oscillation is higher than a certain threshold, the 
system becomes stable. This result is a bit surprising at 
first sight, but can be given an interesting explanation in 
terms of the so-called parametric resonance ||] . 

In this Letter we shall study a system that can be 
viewed as a quantum version of the inverted pendu- 
lum. The system to be considered makes use of down- 
conversion processes interspersed with zones where a lin- 
ear coupling takes place between the down-converted 
photon modes. It is similar to other examples previously 
analyzed in the context of the quantum Zeno effect 
1^ , where the "measurement" is performed by a mode of 
the field on another mode. When the coupling between 
the two modes is large enough, the measurement becomes 
more effective and the dynamics gets stable: this is just 
a manifestation of the quantum Zeno effect, which con- 
sists in the hindrance of the quantum evolution caused 
by measurements. The very method of stabilization of 
the quantum system analyzed here is one of its most 
interesting features and the configuration we discuss is 
experimentally realizable in an optical laboratory. It is 
therefore of interest both for the investigation of the sta- 
ble/unstable borderline for classical and quantum me- 
chanical systems and their links with the quantum Zeno 
effect. 

We consider a laser field (pump) of frequency ujp, prop- 



agating through a nonlinear coupler. The field is con- 
sidered to be classical and the signal and idler modes 
are denoted by a and 5, respectively. We will assume 
that all modes are monochromatic and the amplitudes 
of the fields inside the coupler vary little during an opti- 
cal period (SVEA approximation). The effective (time- 
dependent) Hamiltonian reads (^=1) 

H{t) = UJaO^a + UJbh% + i7int(t), (1) 

where the interaction Hamiltonian is given by 

Hint W - I ^^^t^ ^^^^^ if n < t < Ti + T2 

(2) 

and Hintit + nT) = -ffint(t), with a period T — ti + T2. 
The nonlinear coupling constant F is proportional to the 
second-order nonlinear susceptibility of the medium x(2) 
, fl to the overlap between the two modes ||^ and n — 
0, 1, • • • , is an integer. 

We require the matching conditions cUp — uja + oJb 
and uia = @- The above Hamiltonian describes 
phase-matched down-conversion processes, for nT < t < 
nT + Ti, interspersed with linear interactions between 
signal and idler modes, for nT + ti < t < {n + 1)T. 
Since time is equivalent, within our approximations, to 
propagation length, our system can be thought of as a 
nonlinear crystal cut into A^ pieces, in each of which a, b 
photons are created in a down-conversion process. A sim- 
ilar configuration was considered in Q. Between these 
pieces, no new photons are created by the laser beam, 
but the idler and signal modes (linearly) interact with 
each other, for instance via evanescent waves. See Fig. 0. 
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FIG. I. The system 

By introducing the slowly varying operators a' = 
e"^"*a, h' = e*"''*6, the free part of the Hamiltonian (|l|) 



1 



is transformed away and the Hamiltonian becomes (sup- 
pressing all primes for simplicity) 



H{t) 



Hu = r(at6t +ab) if < i < n, 



(3) 



Hs = 0(a^6 + ab'^) if n < i < ti + T2, 
with H{t+nT) ~ H{t), yielding the equations of motion 



d = —i[a, H], 

In terms of the variables 
1 



i) = -i[b, H]. 



(4) 



x± 



[(a + at)T(6 + fet)], 



p±^~-[{a~a^)Tib^b% 



(5) 



which satisfy the equal-time commutation relations 
= = i, others = 0, the Hamiltonians 

become 



(6) 



They describe two uncoupled oscillators, whose equations 
of motion are 



P± = -«[P±, ^^u] = ±ra;± 

x± = -i[x±,Hs] = ±flp± 
P± = -i[p±, Hs] = T^x± 



x± - r2a;± = 

P± - r^P± = 



x± 

p± - 








(7) 



The first set of equations describes an unstable motion, 
the second set a stable one, around the equilibrium point 
X = p = 0. Notice that the motion of is the 

time-reversed version of that of This is due to 

the fact that the two motions are governed by Hamil- 
tonians with opposite sign in Eq. (|6|). Henceforth, we 
shall concentrate on the variables (a;+,p+) [the stability 
condition for (a;_,p_) is identical]. The solutions are 



P+iTl) 



A, 



x+{0)' 

P+iO) 

cosh(rTi) sinh(rTi) 

sinh(rTi) cosh(rTi) 



(8) 



for the period governed by i?u and 

= A, 



P+{t2) 



fx+iO) 

Sm(flT2) 

sin(ilr2) cos(r2r2) 



^ f C0s(f2T2) 



(9) 



for that governed by Hs . Remember that T = ti + T2 is 
the period of the Hamiltonian H{t) in (^). 

The dynamics engendered by at time t — NT (re- 
member that n = 1, . . . , N) yields therefore 



x+{NT) \ _ ,^/x+(0) 
p+(7VT)J-^ [p+{0) 



A = AsA^. (10) 



These equations of motion have the same structure of a 
classical inverted pendulum with a vertically oscillating 
point of suspension Q, whose classical map is given by 
the product of two matrices Ad = A2A1, with 



Ai 
A2 



cosli(fcir) fcj^ "'^ sinh(fciT 

fcisinh(fciT) cosh(fciT) 

C0s(fc2T) fc2^^sin(fc2T) 

— fc2Sin(fc2T) cos(fc2''") 



(11) 



where the parameters ki and k2 are subject to the phys- 
ical condition ki > k2 > 0. Observe that our system has 
more freedoms: ri and T2 are in general different and the 
parameters fl and F do not have to obey any additional 
constraint. 

The global motion is stable or unstable, according 
to the value of |Tr^| | 2 |]. The stability condition 
|TrA| < 2 reads 

|TrA|/2 = I cos(f2T2) cosh(rri)| < 1. (12) 

This relation is of general validity and holds for any value 
of the parameters n,r and r^. The value of |TrA|/2 is 
shown in Fig. ^). A small-r expansion (the physically 
relevant regime: see final discussion) yields 

1 - (r!V| - FV2)/2 + 0(t4) < 1, (13) 
so that the system is stable for f2r2 > Fri when T2 — > 0. 
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FIG. 2. Stability condition (n2|) in parameter space, a) 
|Tr4|/2 vs Qr2 and Fn; b) Stability (Zeno) region. 

It is interesting to discuss the stability condition just 
obtained for the (x,p) variables in terms of the number 
of down-converted photons. To this end, let us look at 
some limiting cases. [Needless to say, the analysis could 
be done from the outset in terms of Ua and rih and would 
yield an identical stability condition (jl^).] When $1 = 
in (|^) and following equations, only the down-conversion 
process takes place and both — a^a and rif, — b^'b grow 
exponentially with time. There is an exponential energy 
transfer from the pump to the a, b modes. On the other 
hand, if F = and the system is prepared in any initial 
state (except vacuum, whose evolution is trivial) , Ua and 
Ub oscillate in such a way that their sum is conserved (this 
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is due to the property [ria + Ub, Hs] =0). If both and 
r are nonvanishing, these two opposite tendencies (ex- 
ponential photon production and bounded oscillations) 
compete in an interesting way. When Fti > ^T2, in the 
limit Ti —^ 0, the exponential photon production domi- 
nates and there is no way of halting (or even hindering) 
this process: the (external) pump transmits energy to 
the a,b modes. In terms of the {x,p) variables, the sta- 
bility condition (|l3|) cannot be fulfilled and the oscillator 
variables move exponentially away from the origin. The 
opposite situation flT2 > Tti is very interesting and dis- 
plays some quite nontrivial aspects: The motion becomes 
stable and the pump does not transmit energy to the a, b 
modes anymore (the two modes oscillate). 

In general and for arbitrary values of all parameters, 
the action of Hs can be viewed as a sort of measurement 
P, P^|JTl|] , in the following sense: the a mode performs 
an observation on the b mode and vice versa, the pho- 
tonic states get entangled and information on one mode 
is encoded in the state of the other one. For example, 
the condition riT2 — tt/2 yields an "ideal measurement" 
of one mode on the other one, for in such a case the 
states I la, Ob) |0a, lb) evolve into each other. From this 
viewpoint, the stabilization regime just investigated can 
be considered as a quantum Zeno effect in that the 
measurements essentially affect and change the original 
dynamics. In fact, if one considers r2r2 as the "strength" 
of the measurement, by increasing (at fixed flT2) the fre- 
quency of measurements, i.e. by letting ti —^ 0, the 
system moves down along a vertical line in Fig. |^) and 
enters a region of stability (Zeno region) from a region 
of instability. (Notice that it is not necessary to consider 
the Ti — limit ( "continuous measurement" ) in order to 
stabilize the dynamics; there is a threshold, given by the 
curve in Fig. 2 b), at which stability and instability in- 
terchange.) Analogously, at fixed Fti, by moving along 
a horizontal line 17x2 ^ 7r/2 the system enters a region 
of stability because the measurement becomes more "ef- 
fective:" indeed, as emphasized before, flT2 = 7r/2 is a tt- 
pulse condition and leads to a very effective measurement 
of one mode on the other one. It is worth stressing that 
even an instantaneous measurement (projection) can be 
obtained by letting T2 — s- 0, while keeping flT2 finite (the 
so-called impulse approximation in quantum mechanics) , 
and in this case our system yields the standard formula- 
tion of the quantum Zeno effect. 

It is interesting (and convenient from an experimen- 
tal perspective) to consider a single-mode version of the 
Hamiltonian (js]) , in which the down-conversion process is 
replaced by a sub-harmonic generation process (degen- 
erated parametric down conversion). The single- mode 
effective Hamiltonian reads 



^^int — 



H{t)=uja^a + Hint{t), 



(14) 



where the interaction Hamiltonians describing the unsta- 
ble and stable part of the device are 



(r/2)(at2e 
{n/2){a''a + aa^) 



if < t < Ti, 

if Tl < t < Tl + T2, 

(15) 



respectively and i?int(^ + nT) = iJint(i)- By introducing 
the slowly varying operator a' = e*"'"*a, the free part of 
the Hamiltonian ( 14 ) is transformed away and the Hamil- 
tonian becomes (suppressing again all primes) 



i/u = (r/2)(at2 + a2) ifO<i<ri, 

Hs = (f2/2)(a^a -I- aa^) if n < i < n -I- T2, 

(16) 



under which the equation of motion d = —i[a, H] follows. 

In terms of the variables x = {a + a^)/V2,p = —i{a — 
a'^)/\/2 the Hamiltonians read 



-{x^-p"), Hs 



{x'+p^). 



(17) 



These Hamiltonians are identical to the two-mode ver- 
sions ^ describing the decoupled mode (a;_|_,p+), apart 
from the substitution F — F. Hence, the stability con- 
dition is given again by Eq. (^2|), which is even in F. 
Also in this case one can talk of quantum Zeno, but the 
"measurement" is performed by the single mode on itself. 
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FIG. 3. Experimental setup, (a) Possible experimental re- 
alization of the Hamiltonian (|l|)-(|2|). NL, nonlinear crystal; 
Mi,{i = 1,2), semitransparent mirrors; Di, detectors; PBS, 
polarizing beamsplitter, (b) Possible experimental realization 
of the Hamiltonian (|l^)-(|l5|). NLi{i = 1,...,7V), nonlinear 
crystals; PSj, phase shifters; D3, detector. The dotted lines 
indicate which elements are computer controllable. 

It is interesting to discuss a possible experimental re- 
alization of the two situations considered in this Letter. 
The experimental arrangement sketched in Fig. ^(a) cor- 
responds to the two-mode (nondegenerate) case, whereas 
that sketched in Fig. ||(b) to the single-mode (degener- 
ate) case. In Fig. ^a) a type II down-conversion pro- 
cess generates two orthogonally polarized beams of down- 
converted light of the same frequency. The two beams are 
mixed using a polarizing beamsplitter PBS. The stable 
part of the evolution of the system is realized by two 
successive passes of the beams through the beamsplitter. 
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Its reflection coefficient, and hence f2r2, is adjusted by 
rotating it. Mirrors and semitransparent mirrors keep 
sending the beams through the crystal many times. A 
successful stabilization of the unstable system is mani- 
fested in the decrease of the rate of photon registrations 
at detectors Di, D2 at a certain position of the beam- 
splitter PBS. A different setup is sketched in Fig. ^(b), 
where N processes of subharmonic generation take place 
in N nonlinear crystals with controlled phase shifters 
in between them. For appropriately chosen phase shifts 
9i = {^T2 + Ci) mod 27r, where Ci are — 1 phase shifts 
intrinsic to the actual experimental arrangement (given 
by distances between crystals, etc.), the generation of the 
subharmonic wave is suppressed. 

In order to give a reasonable estimate of the value of 
the coupling constant T, consider that, due to the cor- 
respondence principle, the gain of classical and quantum 
parametric amplifiers must be the same; therefore one 
can use the well-known classical formula for the nonlin- 
ear coupling parameter Fc governing the space evolution 
inside the nonlinear medium, which in MKS units reads 

= \x{'if^a^bip. (18) 

Here 77 is the impedance of the medium, x(2) is the 
second-order susceptibility, and are the frequen- 
cies of modes a and 6, respectively, and Ip is the inten- 
sity of the pump beam. The following numerical values 
could be typical for a performed experiment: rj « 2200, 
X(2) w 2 X 10-23 CV-2, = W6 w 3 X IQi^s-i and 
Ip ~ 10''^Wm^2. Hence the nonlinear coupling parame- 
ter is of the order of Fc « O.lm"^. Reasonable lengths of 
nonlinear crystals are of the order of I « lO^^m, so that 
the dimcnsionless product of interest can be estimated to 
be about 

Fti = FcZ « 0.001. (19) 

This means that the down-converted beam(s) ought to 
pass the nonlinear region many times in order to show 
an explosive increase of its (their) intensity (ies). This 
could be achieved by placing the nonlinear crystal in a 
resonator as shown in Fig. ^(a). However, in order to 
observe a significant change of the dynamics of the pro- 
cess in question due to the performed stabilization, a few 
passes might already turn out to be sufficient. 

In conclusion, we have discussed a striking quantum- 
optical analogue of a well-known classical unstable sys- 
tem. By interspersing the nonlinear regions with regions 
of suitably chosen linear evolution, the global dynamics 
of our system can become stable and the generation of 
down-converted light can be strongly suppressed. This 
behavior has an interesting interpretation in terms of the 
quantum Zeno effect: by increasing the "strength" of the 
observation performed by the a mode on the b mode and 
vice versa, in the sense discussed before, the evolution is 



frozen and the system tends to remain in its initial state. 
This phenomenon is somewhat counterintuitive: in the 
setups in Fig. ||, even though the beams are forced to go 
through the crystal many times, no exponential photon 
production takes place. The experiment seems feasible 
and its realization would illustrate an interesting aspect 
related to the stabilization of a seemingly explosive be- 
havior. 
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